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Abstract 

In this paper fractional generalization of Liouville equation is considered. We derive 
fractional analog of normalization condition for distribution function. Fractional gener- 
alization of the Liouvile equation for dissipative and Hamiltonian systems was derived 
from the fractional normalization condition. This condition is considered considered 
as a normalization condition for systems in fractional phase space. The interpretation 
of the fractional space is discussed. 
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We call a fractional equation a differential equation that uses fractional deriva- 
tives or integrals. Fractional derivatives and integrals have found many appli- 
cations in recent studies of scaling phenomena. We formulate fractional analog 
of main integro-differential equation to describe some scaling process - Liouville 
equation. Usually used for scaling phenomena Fokker-Planck equation can be 
derived from Liouville equation. Therefore it is interesting to consider a frac- 
tional generalization of the Liouville equation. To derive fractional equation for 
a distribution function we must consider a fractional analog of the normalization 
condition for distribution function. Most of the fractional equations for distri- 
bution function does not use correspondent normalization condition. Therefore 
these equations (with fractional coordinate derivatives) can be incorrect equa- 
tions. In this paper fractional Liouvile equation for dissipative systems is derived 
from the normalization condition. The coordinate fractional integration for this 
normalization condition is used. 
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1 Introduction 



Fractional derivatives and integrals have found many applications in recent studies of scaling 
phenomena PjJ El HI El • The main aim of most of these papers is to formulate fractional 
integro-differential equations to describe some scaling process. Modifications of equations 
governing physical processes such as the Langevin equation jH], diffusion equations, and 
Fokker-Plank equation have been suggested [7j-^3| which incorporate fractional derivatives 
with respect to time. It was shown in Ref. [21 that the chaotic Hamiltonian dynamics of par- 
ticles can be described by using fractional generalization of the Fokker-Planck-Kolmogorov 
equation. In Ref. [Hj, coordinate fractional derivatives in the Fokker-Planck equation were 
used. 

It is known that Fokker-Planck equation can be derived from Liouville equation ^Hl EE] • 
Therefore it is interesting to consider a fractional generalization of Liouville equation and 
Bogoliubov hierarchy equations. We call a fractional equation a differential equation that 
uses fractional derivatives or integrals. 

To derive fractional equations for a distribution function we must consider a fractional 
analog of the normalization condition for distribution function. Fractional Liouville equa- 
tion for dissipative systems is derived from the normalization condition. In this paper, the 
coordinate fractional integration for normalization condition is used. This condition is con- 
sidered considered as a normalization condition for systems in fractional phase space. If 
any fractional equation for distribution function does not use correspondent normalization 
condition, then this equation (with fractional coordinate derivatives) can be incorrect. 

In Sec. 2 the normalization condition for distribution function and notations are consid- 
ered. In Sec. 3 we derive the Liouvile equation from the normalization condition. In Sec. 
4 the physical interpretation of fractional normalization condition is considered. Finally, a 
short conclusion is given in Sec. 5. 

2 Normalization condition 

Let us consider a distribution function p(x,t) for x G R 1 . Let p(x,t) G L 1 (_R 1 ), where t is a 
parameter. Normalization condition has the form 



where y G (— oo, +oo). 

Let p(x,t) G LplR 1 ), where 1 < p < 1/a. Fractional integrations on (— oo,y) and 
(y, +oo) are defined [17J by 




p(x, t)dx = 1. 



This condition can be rewritten in the form 






(2) 
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Using these notations, Eq. has the form 

Vlp)(y,t) + (llp){y,t) = l. 
Using definitions (j2J) and (J3J) we can get the fractional analog of normalization condition 

(I«p)(y,t) + (I«p)(y,t) = l. 
Equations (j2J) and (jSJ) can be rewritten in the form 

(/» (y, «) = tvT / * a ~ V(y T t)dx. (4) 
1 a Jo 



This leads to the normalization condition 

x a ~ l [p(y — x, t) + p(y + x, t)\dx = 1. (5) 



l 

1 ' _a-lr 







r(a) 

If we denote 

p(x,t) =p(y-x,t) +p(y + x,t) (6) 

and 

r(a) 

then condition (jlj has the form 

/■oo 

/ p(x,t)dp a (x) = 1. (8) 
Jo 

Note that substituting y = rf in (JBJ, we get the sum 

p(x t , t) = p(ct - x t , t) + p(ct + x t , t). 
This sum can be considered as a sum of right and back waves of the distribution functions. 

3 Liouville equation 

Let us consider a domain B for the time t = 0. In the Hamilton picture we have 



p(xt,t)dp, a (x t ) = p(x ,0)dp a (x ). 

B t J Bo 

Using the replacement of variables Xt = Xt(xo), where xq is a Lagrangian variable, we get 



p(x t ,t)xf ——dx = p(x ,0)xq dx . 

B OX JBo 

Since B is an arbitrary domain we have 

p(x t ,t)dp a (x t ) = p(x ,0)dp a (xo), 



or 

p(x t ,t)xr 1 ^ = p(xo,OK- 1 . 
Differentiating this equation in time t, we obtain 

dp(x t ,t) a _ x dx t t f^^dxt 

-X-t 



dt dxo dt \ dxo ■ 

or 



dt 

where d/dt is a total time derivative 



(1/,(r ''Kn a (x t ,t)p(x t ,t) = o, (9) 



d d „ d 

— = h F t . 

dt dt dx t 



The function 

describes the velocity of (phase space) volume change. Equation © is a fractional Liouville 
equation in the Hamilton picture. If the equation of motion has the form 

then the function Q Q is defined by 
As the result we have 

The normalization in the phase space is derived by analogy with a normalization in the 
configuration space. The fractional normalization condition in the phase space 



o JO 



p(q,p,t)dp a (q,p) = 1, (10) 



where dp Q (q,p) has the form 



dq a A dp a (qp) ' 1 
dp a {q,p) = dp a {q) A d/i a (q) = = dq A dp. (11) 

(al ya)) 1 yet) 

The distribution function p(q, p, t) in the phase space is defined by 
p(q, P, t) = p(q' -q,p' ~ p, t) + p(q' + q,p' - p, t) + p(q' -q,p' + p, t) + p{q' + q,p' + p, t) . 
Let us use the well known transformation 

dq t A dp t = {qt,Pt}odq A dp , (12) 
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where {qt,Pt}o is Jacobian which is defined by the determinant 

I„ „\ - dct d{ - quPt) - det ( d ^t/9qm dq kt /d Pm \ 
{qUPt} ° ~ %o,Po) ~ { d Pkt /dq m d Pkt /d Pl0 ) ■ 

Using ptd/i a (q t ,Pt) = podp, a (q ,Po), we get the relation 

Pt „ 2 , r dq t A dp t = p dq A dp . (13) 

Using (|T2|) . we have condition ()13|) in the form 

PMtPtT~ l {quPt}o = {qoPoT^Po (14) 
Let us write condition ()14|) in more simple form 

Pt{q?,pno = a 2 (q p o r- 1 po. (15) 
The time derivatives of this equation lead to the fractional Liouville equation 

dp(q t ,p t ,t) 



dt 

where the function Q Q is defined by 



+ n a (qt,p t ,t)p(q t ,Pt,t) = 0, (16) 



n a (q t , Pt ,t) = {gf,j£}o 1 ^{tf,rf}o = j t ln{q?,pt} Q . (17) 



In the usual notations we have 



tt*(q t ,Pt,t) = ^-ln det™' 1 ' '''' 



dt ' d(q ,Po) ' 

Using well-known relation In det A = Sp In A, we can write the a-omega function in the 
form 

1 dt + dt la ' 

where 

dAdB _ dA dB 

In the general case (a ^ 1) the function f2 a is not equal to zero (Q a ^ 0) for Hamiltonian 
systems. If a = 1, we have Q a ^ only for non-Hamiltonian systems. 
It is easy to see that any system which is defined by the equations 



dq" p t dpf 



dt m dt 



filt), (19) 



has the a-omega function equal to zero Q a = 0. This system can be called a fractional 
nondissipative system. For example, a fractional oscillator is defined by the equation 



The a-omega function can be rewritten in the form 



_ x dq t _ x dp t \ ( dq t dp t 



Q a (q t , Pt ,t) = (a- l)[q-^+ p -^j + {~^,Pt}x + {ft, -£h, (21) 
where {., .}i is usual Poisson bracket. If the Hamilton equations have the form 

^ = G(q t ,p t ), ^ = F(q t ,p t ), (22) 

then the a-omega function is defined by 

Q a {q,p) = (a- l)(q- 1 G(q,p)+p- 1 F( q ,p)) + {G,p}i + {q, F} t . (23) 



This relation allows to derive fl a for all dynamical systems (|22j) . It is easy to see that the 
usual nondissipative system 

dq t Pt dp t 

-77 = — , -7r = f{Qt), 24 
dt m dt 

has the a-omega function 

M a (q,p) = (a - l)(mqpY 1 (p 2 + mqf(q)) 

and can be called a fractional dissipative system. For example, the linear harmonic oscillator 
(/(<?) = -muj 2 q) 

&a(q,p) = (a — l)(mqp)^ 1 (p 2 — m 2 u> 2 q 2 ), 
is a fractional dissipative system. 

4 Interpretation 

The fractional normalization condition can be considered as a normalization condition for 
the distribution function in a fractional phase space. In order to use this interpretation we 
must define a fractional phase space. 

The first interpretation of the fractional phase space is connected with fractional dimen- 
sion. This interpretation follows from the well-known formulas for dimensional regulariza- 
tions 

27r n / 2 r°° 

p(xKx = TW2) I p{x)xn ~ dx - (25) 

Using Eq. (|25j) . we get that the fractional normalization condition (JHJ) can be considered as 
a normalization condition in the fractional dimension space 



r(q/2) 
2W 2 ri 



a 



p(x,t)d a x=l (26) 



up to the numerical factor r(a/2)/(27r Q / 2 r(a)). 

The second interpretation is connected with the fractional measure of phase volume. The 
parameter a defines the space with the fractional phase volume 



V a = d[i a {q,p). 
Jb 
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It is easy to prove that the velocity of the fractional phase volume change is defined by 



n a (q,p,t)dn a (q,p). 



dt Jb 

Note that the volume element of fractional phase space can be realized by fractional exterior 
derivatives [T5] 

n Qo n Qa 

d a = V dql— -— + V dpi— — — , (27) 

in the form 

dq a A dp- = (t^T^) + f^ha)y liqPr ~ ldaq A d> (28) 

The system can be called a fractional dissipative system if a fractional phase volume 
changes, i.e., Q a ^ 0. The system which is a nondissipative system in the usual phase space, 
can be a dissipative system in the fractional phase space. The fractional analog of the usual 
conservative Hamiltonian nondissipative system is defined by the equations 

M = *M, M = A(9) . (29) 

dt m dt 

The usual nondissipative systems (J24j) are dissipative in the fractional phase space. 

In the general case, the fractional system is a system in the fractional phase space. We 
shall say that a system is called a fractional system if this system can be described by the 
fractional powers of coordinates and momenta, 

„a |„ \a I la 

Qk - \Qk\ , P k - \Pk\ , 

where k — 1, n. 

The fractional systems allow to consider the interpretation of the fractional normaliza- 
tion condititon which is used to derive the fractional Liouville equation. The fractional 
normalization condition for the distribution function can be considered as a normalization 
condition for the systems in the fractional phase space. 

The Hamilton equations for the fractional system have the form 

d -§ = P l d § = F k {q^p«). (30) 
at m dt 

Obviously, that the equation dp k /dt = F k can be rewritten in the fractional form. Mul- 
tiplying both sides of this equation by ap a ~ l , we obtain dp^/dt = ap a ~ 1 F k . However the 
equation dq k /dt = p k /m cannot be rewritten in the fractional form dq^/dt — p%/m. 
The fractional conservative Hamiltonian system is described by the equation 

dqt = dH(q a ,p a ) dp% dH(q a ,p a ) 

dt dpi ' dt dqt ' 1 ' 

or 



^ = {fi,H} a , ^§ = {P a k,H}a- (32) 



7 



Here H = H(q a ,p a ) is a fractional analog of the Hamiltonian function 



n 2a 

H(q«,p") = J2^ + U(q a ). (33) 



The fractional system can be considered as a nonlinear system with 



n i 

H(q a ,p a ) = £ -g kl (q,p)PkPi + U(q). (34) 
fc=i ^ 



Note that the Hamiltonian (}3*4"j) defines a nonlinear one-dimensional sigma-model 
with metric 

It is easy to see that fractional systems (J33|) can lead to the non-Gaussian statistics. 
The interest in and relevance of fractional kinetic equations is a natural consequence of the 
realization of the importance of non-Gaussian statistics of many dynamical systems. There 
is already a substantial literature studying such equations in one or more space dimensions. 

Note that the classical (nonlinear) dissipative systems can have canonical Gibbs distri- 
bution as a solution of stationary Liouville equations for this dissipative system [221 • Using 
the methods (221, it is easy to prove that some of fractional dissipative systems can have 
fractional canonical Gibbs distribution 

p(q,p) = Z(T)exp — ' J ' 



kT 

solution of the fractional Liouville equations 



i^ + M^ M ^°- (35) 

Here the function H(q a ,p a ) is defined by (|33|). 



5 Conclusion 

Derivatives and integrals of fractional order have found many applications in studies of scal- 
ing phenomena PQ |2l El HI El • In this paper we formulate fractional analog of main integro- 
differential equation to describe some scaling process - Liouville equation. We consider the 
fractional analog of the normalization condition for the distribution function. Fractional Li- 
ouvile equation for dissipative systems is derived from the normalization condition. In this 
paper, the coordinate fractional integration for the normalization condition is used. This 
fractional normalization condition can be considered as a simulating unconventional envi- 
ronment for systems with fractional dimensional phase space or phase space with fractional 
powers of coordinates and momenta. Note that the adoption of fractal formalism yields 
properties that the ordinary formalism would produce only in the case where the system is 
made non-Hamiltonian by the presence of an environment, whose influence can be mimicked 
by means of friction, for instance. 
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Suggested fractional Liouville equation allows to formulate the fractional equation for 
quantum dissipative systems |2Hj by methods suggested in Refs. j2H 123 • In general, we 
can consider this dissipative quantum systems as quantum computer with mixed states ■ 
These dissipative quantum systems can have stationary states [27] ■ Stationary states of dis- 
sipative quantum systems can coincide with stationary states of Hamiltonian systems [2H] ■ 

This work was partially supported by the RFBR Grant No. 02-02-16444 
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